Nucleation and growth of second phase precipitates under non-isothermal conditions * 
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The Langer-Schwartz equations for precipitation are formulated to calculate nucleation, growth 
\ and coarsening of second phase precipitates under non-isothermal situations. A field-theoretic 

steady-state nucleation rate model is used in the analysis. The field-theoretic nucleation rate is 
compared with the classical nucleation rate relation derived from the Fokker-Planck equation. This 
^ \ integrated model is used to simulate a bcc-to-hcp phase quenching and subsequent annealing in 

O , the hep phase of a dilute zirconium alloy, where the mean precipitate size, number density and the 

' degree of supersaturation are calculated as a function of time. The influence of cooling rate on the 

aforementioned parameters is evaluated. A lower cooling rate results in larger precipitates with a 
00 , smaller number density in concordance with observations. 

• 1— I \ 

^ '. I. INTRODUCTION 

C/2 , 

?H ■ Various important properties of engineering alloys, such as their mechanical strength and toughness, creep and 

SI corrosion resistance, magnetic and superconducting behaviors, are basically governed by the presence and character- 
. . istics of second phase precipitated particles (SPPs) In many material processing methods, the SPP characteristic 
is controlled by quenching the material from a high temperature phase, where the constituents of the SPPs are in 
soHd solution, to a low temperature phase, in which precipitation occurs and a distribution of SPPs is observed. This 
^ • process embraces such phenomena as nucleation, spinodal decomposition, late stage growth and coarsening, which 

' are broadly described by the framework of the kinetics of first-order transitions 0. 
P5 ■ For example, when a binary alloy system is quenched from an equilibrium state (phase) to a non-equilibrium state 
O I inside a coexistence curve of its phase diagram, the quenched system gradually transforms from the non-equilibrium 
, state to an equilibrium state, consisting of two coexisting phases. This transformation occurs by time-dependent 

spatial fiuctuations of the local concentration of one of the two component species. Commonly, two different kinds of 
I . instabilities are considered; in one case the soHd solution experiences a shallow quench in the metastable region, in 

■ another, it is quenched deeply into the unstable region of the miscibility gap [ij. In the former case, the instability is 
due to the formation of stable nuclei via localized composition fiuctuations, requiring to surmount an energy barrier 
that is characterized by an incubation period. This corresponds to a metastable phase and the transformation is 
by nucleation and growth. In the latter case, the instability is caused by the formation of non-localized, small- 
amplitude, spatialhy extended composition fiuctuations that grow spontaneously in amplitude as the time after the 
quench proceeds For a binary system, this phenomenon is called spinodal decomposition, for which the quenched 

. state Hes within the boundary between unstable and metastable states, the so-called spinodal curve. 
' Langer and Schwartz 0| have developed a detailed theory describing the nucleation, growth and coarsening of 
[ droplets in metastable, near-critical fiuids. Wendth and Haasen j3| adapted this theory for solid solutions and used 
2 ■ it to interpret the precipitation of 7'-Ni3Al SPPs in Ni-14 at.%Al. Later, Kampmann and Wagner introduced a 
non-linearized Gibbs-Thomson relation into the Langer-Schwartz equations in order to account for cases where the 
surface tension of precipitates and the supersaturation are large. They also included an empirical time-dependent 
Q ' factor on the steady-state nucleation rate to include transient nucleation kinetics that has been observed in many 
Q . experiments. Kampmann and Wagner used the model to interpret and explain observations on /?'-Cu4Ti particles in 
' Cu-1.9 at.%Ti alloy. 

. ^ ' In the present paper, we consider homogeneous nucleation, growth and coarsening of SPPs under non-isothermal 

■ conditions. The Langer-Schwartz equations for the precipitation kinetics are formulated in a way that can be used 
H [ to calculate nucleation and growth of SPPs under time-varying temperature situations such as quenching. A field- 
. . . theoretic steady-state nucleation rate model is used in our evaluation. Moreover, the basic models for growth and 

coarsening are presented, followed by a computation of an isothermal experiment on a copper-cobalt alloy phase 
transformation, and a simulation of a body-centered cubic (bcc) to hexagonal close-packed (hep) phase quenching and 
subsequent annealing in the hep phase of a dilute zirconium alloy. 
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* This is an expanded version of the paper presented in the International Conference on Sohd-Solid Phase Transformation in Inorganic 
Materials 2005 (PTM 2005), May 29 - June 3, 2005, Phoenix, Arizona, USA. 
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II. MODELS 



The models considered here describe the phenomena of homogeneous nucleation, growth and coarsening of second 
phase precipitates in non-isothermal conditions in an integrated fashion. The precipitates are assumed to be spherical 
and at each instant are in local equilibrium with their surrounding; therefore the concentration near the SPP / matrix 
interface is determined by the Gibbs-Thomson boundary condition. A key variable defining the state of the system 
is the supersaturation; here defined as a; = ln(C/Coo), where C is the solute concentration in the matrix, which is 
time-dependent, and Coo the solubility limit or the equilibrium solute concentration, which is strongly temperature 
dependent. The Langer-Schwartz theory of precipitation basically comprises a nucleation model and a growth model. 
In this section we outline these models and describe our formulation of the theory for non-isothermal applications. 



A. Nucleation 

1. Classical Theory 

Consider that a spherical (or circular) particle with radius R is emerged from a metastable state. The particle free 
energy consists of the surface and the bulk energy parts, 

F(i?) = 5rf(ai?^'-i-^i?''), (1) 

where d is the spatial dimensionality, 5*3 = 47r and 5*2 = 27r, ct the surface energy, and /ie the difference in free energy 
density (per unit volume) between the metastable and stable phases. The function F{K) goes through a maximum 
at a critical radius R = Rc = {d — l)a/^e- The height of the nucleation energy barrier is given by: 

Fc ^ F{Rc) = ^aRi-\ (2) 

The probability of formation of a particle of radius Ris P ^ exp [—I3F{R)], where P = l/fc^T with ks Boltzmann's 
constant and T the temperature. Thus the critical nucleus, which maximizes F, is the least probable one. Once the 
particle forms, it will grow to reduce its free energy. 

The basic ingredient of the considered model is the nucleation rate of the secondary phase precipitates. In the 
standard formulation, attributed to Ya. B. Zel'dovich 0, the nucleation process is described by a Fokker-Planck 
equation: 

dn d f^dn\ dF\ 

where n — n{R,t) is the number of SPP in the radius interval [R,R + dR] per unit volume at time t, B ~ B{R) a 
nucleation size diffusion coefficient and F the SPP free energy. After a transient time to, a constant nucleation rate of 
SPPs with radii larger than Rc prevails, J = const. More concretely, if R » Rc, n{R,t)dR = Jdt is the number of 
growing newly emerging particles in a time interval dt per unit volume and is independent of i? . The steady-state 
solution of Eq. lO, n = ns{R) is obtained by integration, viz., 

POO 

n,(i?) = Je-™-) / dri?(r)-ie'5^W, (4) 
Jr 

where the condition ns{R) ^ as i? ^ oo was imposed. Next, noting that the steady-state solution of the Fokker- 
Planck equation is determined by the Boltzmann distribution, i.e., Us = n^exp [—f3F{R)], an expression for J can be 
obtained, when R is near Rc, giving (AppendixEl 

J = (2^)-l/2(d - 1)1/2 A.(5d/3fT)-l/2i?-('^+l)/2n, 

where I?„ = SdPcrRf^^Bc is the critical nucleation diffusivity with Be = B{Rc). Relating the critical radius to the 
supersaturation x and the capillary length £ according to: Rc — i/x, where £ = 2(3ava and Va the atomic volume of 
the precipitate, Eq. ^ is expressed in terms of the supersaturation, i.e.. 
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where 

Relation ((^ is the classical expression for the nucleation rate formulated in d spatial dimensions as a function of the 
supersaturatio n. The scalin g in the critical region suggests that the Boltzmann prefactor is ~ Whence, 
setting A = \/ {d — l)/27r(i, Eq. © can be written as 



2. Field Theoretic Approach 

The nucleation rate in the field-theoretic formulation Q begins from the relation that links the steady-state nucle- 
ation rate with the imaginary part of the free energy 

J=^^[F{^e)\. (9) 

where k is a kinetic factor related to the diffusivity via k — {d—l)Dni/Rc- Relation jHJ has been evaluated by various 
investigators and more carefully by Giinther, Nicole and Wallace which for d = 3 takes the form (Appendix El 

T~i a / \ 2/3 



where B — 6^/(2887r\/3). Equation ifTTHl is strictly valid for low supersaturations, i.e., when x/xq « 1. Langer and 
Schwartz [4} have heuristically extended this expression for large values of x according to 

_ -LJ fi X f) I X \ I X 



The behavior of the aforementioned models for the nucleation rate, Eqs. ijHJ, and are illustrated in Fig. 
2] for d = 3, where the scaled nucleation rate K is plotted against the scaled supersaturation x/xq (Eq. ^ is scaled 
by ADnXo/i'^+^, while Eqs. CHI and ^ are scaled by BDnX^/£^). 



B. Growth and Coarsening 

1. Theory 

The details of the Langer and Schwartz theory have been clearly described in their original work and reviewed 
in a number of publications (Q, 0|, We only summarize its basic ingredients in order to fix notation. The 

basic equations of the Langer-Schwartz theory are those of the Lifshitz and Slyozov [JJJ , supplemented by a source 
of droplets (spherical precipitates) which describes the nucleation event. The continuity equation for the droplet 
distribution function n — n{R,t) is 



dn d 



dR 
dt 



where j{R) is the distributed nucleation rate of precipitates, J — j{R)dR, and as before J stands for the nucleation 
rate per unit volume. Langer and Schwartz made a simplifying assumption that only precipitates with R > Rc are 
to be counted as part of the second phase. The SPP number density N is given by iV = n{R)dR with their mean 

radius defined as R — n{R)RdR. 

The growth rate of the spherical nucleus is controlled by the rate the solute atoms supplied to the precipitate / matrix 
interface via diffusion [lal 



dR D f C{t) - Cr 
dt R \ dp — Cfl 



(13) 
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FIG. 1: The scaled nucleation rate vs. scaled supersaturation for d = 3, comparing predictions of Eq. the classical model, 
Eq. linil . the GNW model, and Eq. the LS model. The scaling factors used are defined in the text 



Here, D is the diffusivity of the solute in the matrix, C the matrix average solute concentration, Cp the solute concen- 
tration within the precipitate, assumed to be uniform, and Cr the local solute concentration at the precipitate/matrix 
interface. The latter concentration is related to the solubility limit Coo, a temperature dependent quantity, and the 
curvature of the precipitate/matrix interface through the Gibbs-Thomson relation Cr = Cqo exp Also, the 

equation for conservation of matter must be satisfied, viz.. 

For non-isothermal conditions, where temperature is varying with time, all the temperature dependent parameters 
appearing in the foregoing equations must be considered and treated properly. The resulting Langer-Schwartz dif- 
ferential equations for N and R, applicable to both isothermal and non-isothermal conditions, can be expressed by 



dN _ ffi - gi/2 + 33/1 . dR _ ./i - /1.92 + /sffi 
dt 1 - /2 - 3253/3 ' dt 1 - .92 - /2/3g3 ' 

where the functions gi, 52, gs and /i, /2, fa are 



(15) 



bm (I IdlnCooN dT ^^^^ 



X( 



{R-Rc) \T X dT J dt 



-NR -N^ 

D ( eHJ^^^_C^\ J , be ^1 din Coo 1\ dT 

h-NT, f.^^T, T^ Co) ^^^^ 



3 



In the above expressions, SRc is the width of the size distribution and 6 = 0.317014 for R < 3Rc/2 and 6 = otherwise 
01 • In computations, we assume SRc = at with a, a positive constant a < 1, taken as a free parameter. It is noted 
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that the particular value of b set by Langer and Schwartz reproduces the familiar Lifshitz-Slyozov coarsening law 
as t ^ oo 0, 01 . The initial conditions for the equations in (QSl are 

N{t = to) = No, (20) 
— £ 

R{t = to) = Rc + SRc = - + ae, (21) 

where the starting time to is determined from No = /J ° J{T)dT; and No defines the lowest particle density of practical 
interest, here treated as a model parameter. 

2. Numerical Method 

The nonlinear ordinary differential equations lfT5)l need to be evaluated numerically. In evaluations of the right- 
hand-side expressions of equations in lfT5|l , we make use of correlations for the solubility limit Coo and the diffusivity 
D, which are temperature-dependent material properties. We have used the MATLAB programming environment to 
solve these differential equations; more specifically the MATLAB solver 0DE15S ^3|- This is a variable order solver, 
which is intended to solve stiff systems of ordinary differential equations. The solver is invoked at each time step of 
the time-temperature history, and the evolution of N and R during the time step is computed with 0DE15S for each 
mesh point, based on the local temperature at beginning and end of the time step. In the solution procedure, the 
temperature is assumed to vary linearly with time during the time step. 

III. APPLICATION 
A. Isothermal experiment 

The methods outlined in the foregoing section have been verified against a number of measurements made 
in isothermal conditions for binary precipitates such 7'-Ni3Al in Ni-14 at.%Al and Co in Cu-Co alloys 0, ITR 
H^ . Here, we only report the results of our computations for cobalt precipitates in Cu-2.7 at.%Co alloy, for which 
isothermal annealing tests (at 823 K) were carried out by Wendth and Haasen The considered experiment 

involved measurements of the SPP number density as a function of annealing time and also the time variation of 
precipitate mean size was determined. These properties were measured by transmission electron microscopy (TEM) 
and atom probe field ion microscopy. 

In order to simulate this test, we used the solubility limit for cobalt in copper proposed by Servi and Turnbull 0| 
Coc = 712.85 X 10~^*^^/^, where T is the absolute temperature and Coo is in atomic percent. For the diffusivity of 
cobalt in copper, we used the correlation by Dohl et al. [20j, D = 4.3 x 10~^e~^^^'^®/^, where D is in m^s~^. The 
matrix/SPP interface energy a was in our analyses set to 0.22 Jm~^, following the evaluation of data by Stowell [2ll |. 
The second phase particles in these alloys are of pure cobalt, which means that Cp is 100 at%. Both the particles and 
the matrix have a face-centered cubic (fee) crystal structure. Moreover, the molar volume Vm of cobalt is 6.7 x 10"^ 
m'^mole"^. In analyses, we have set the model parameter No in Eq. H2()ll to 1 x 10^° m~^ and a in Eq. Ij21|l to unity. 

The calculated evolution of SPP number densities are compared with data in Fig. [2l The results are shown for 
two cases: For the nominal case, the nucleation rate J, given by Eq. lfTT|l . has been used. For the scaled case, the 
coefficient B in Eq. (Illll has been scaled by a factor 0.05, in order to obtain a best fit to the data. For this test at 
823 K, shown in Fig. |2l the measured SPP density reaches a peak at about 300 seconds, after which coarsening starts 
to dominate the picture. At this stage, large particles will grow at the expense of smaller ones, as a consequence 
of the Gibbs-Thomson relation, see section IIIB II The smallest particles will ultimately dissolve, which makes the 
number density decrease with a rate proportional to ~ 1/t. Our model captures this behavior fairly well, but the 
calculated peak in the SPP density is reached too early, in comparison with data. For Cu-1.0 at%Co alloy, LeGoues 
and Aaronson ||1^] reported a delay time of about 130 s in their study, whereas Haasen and Filler ^] reported that 
the delay time was no more than 30 s in the study made by Al-Kasab. Since we have utilized a steady-state theory 
and assumed a homogeneous nucleation rate, we could have overlooked the effect of the delay time, td, and hence get 
a discrepancy between the observed temporal data and calculated values; cf. Appendix [0 for a rough estimation of 
this effect. 
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FIG. 2: Calculated number density of Co second phase particles in Cu-2.7at%Co at 823 K, in comparison with data from [T^ . 
The calculated nucleation rate, given by Eq. lllll . is scaled by a factor 0.05 in order to fit the data. 



B. Non-isothermal simulation 



In this subsection, we present the results of a simulation of a heat treatment made on a cylindrical specimen. The 
cylinder is made of Zircaloy-2 (Zr-l.4Sn-0.12Fe-0.09Cr by wt.%) and has a diameter of 25 mm. The heat treatment 
consists of quenching the sample from a high temperature (1323 K) body-centered cubic /? phase in water, which is 
kept at room temperature. After the quenching, the cyHnder is subjected to two subsequent annealing steps in Ar gas 
at 838 K (hexagonal-closed packed a phase) for durations of 1.0 h and 1.5 h, respectively (Fig. EJ. The prominent 
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FIG. 3: Temperature history for the annealing steps; the cooling and heating rates for the steps are ±5.4 Ks 



second phase precipitate in Zircaloys is Zr(Fe,Cr)2 (Laves phase with hexagonal structure) formed at around 1200 
K [23|- Here, we calculate the nucleation and growth of this kind of precipitates under the aforementioned heat 
treatment (quenching plus anneahng) using the models described in the foregoing section, which are implemented in 
a computer program [231 • In particular, we have selected the nucleation rate equation llll|l in our evaluation. Since 
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the diflusivitv of Fe and Cr are similar and solubility data (temperature vs. concentration) are available for (Fe+Cr) 
in Zircaloy [23|, we have treated the precipitate as a binary compound Zr and (Fe+Cr) with initial solute (Fe+Cr) 
concentration of Co = 2010 wppm (weight parts per million). The SPP molar volume is Vm — 9 x 10~^ m'^mole^^ 
and its composition Cp = 5.4 x 10^ wppm. We have assumed that the diffusivity during nucleation is identical to 
that during growth. The input model parameters in our evaluations are as follows. We set the effective diffusivity for 
(Fe+Cr) in Zircaloy as D{T) ~ 1.473 x io~6g-i5930/T ^ -^jigj-g x is the absolute temperature and D is in m^s^^, and 
the surface tension of the precipitate cr — 0.25 Jm~^. The lower cut-off limit for the particle density, No in Eq. (12(1)1 . 
is set to 1 X 10^° m~^, and the model parameter a in Eq. Ij21|l to 0.25. The motivations for selecting these values are 
discussed in [l5| . 

The results of our non-isothermal evaluation are presented in a number of diagrams in Fig. 31 which shows (i) the 
mean precipitate radius, (ii) the precipitate number density, (iii) nucleation rate, (iv) the matrix supersaturation, (v) 
the matrix solute concentration (Fe+Cr), and (vi) temperature evolution in the first ten seconds of the heat treatment. 
The computation output corresponds to the center of the cylinder and 0.1 mm from the outer surface. The cooling 
rates at these locations, during quenching, can be noted from the temperature vs. time diagram in Fig. 31 

C. Discussion 

From the results displayed in Fig. 31 we note that there is a significant difference in calculated precipitate radius 
at the different positions in the cylinder, directly upon {3 quenching, but this difference gradually diminishes under 
subsequent annealing. The initial difference in the radius of the precipitates located in the cylinder central region and 
those close to the outer surface follows from the difference in the cooling rate, as is seen in the temperature diagram in 
Fig. ^ According to our computations, the slow cooling at the cylinder central part results in appreciable precipitate 
growth already under the quenching phase. This behavior is not predicted at the billet outer surface. 

The calculated time variation of the precipitate number density shown in Fig. 31 displays a remarkable difference 
between the central and periphery of the cylinder directly on quenching. The calculated number density at the outer 
surface is initially about several orders of magnitude larger than that at the central part, but this large difference grad- 
ually disappears under subsequent annealing as a result of coarsening. Figure 31 also shows the calculated nucleation 
rate as a function of time, evaluated at the considered locations in the cylinder. As can be seen, the nucleation takes 
place under a very short time span, as the temperature drops below 1118 K and the matrix becomes supersaturated 
(C > Coo). The nucleation rate pulse width, as can be seen in the diagram, is a couple of ms at the billet outer 
surface, whereas the corresponding time span is much larger in the central part, i.e., around 50 ms. The calculated 
matrix supersaturation, x = ln(C/Coo), is also shown as a function of time in the figure. The spikes seen in the 
diagram for x correspond to the temperature drops during the annealing cycles (Fig. , which reduce Coo • Note also 
the sharp dive in the matrix solute concentration due the temperature dip. The results of the computations presented 
here, i.e., the precipitate radius and number density, are in qualitative agreement with the observations made on a 
similar kind of heat treatment of this material [23 |. 

As alluded in section ITlI Al we have supposed a steady-state nucleation model. Such a treatment does not provide 
information regarding instantaneous SPP size distribution nor on the nucleation rate prior to reaching steady state. 
Some workers postulate a nucleation relation in the form J{t) = J (oo) exp {—td/t) to fit metallurgical data P|. 
However, this typ e of approach is simplistic and a more rigorous treatment to solve the time-dependent Fokker-Plank 
equation (pj or extend the field theoretical approach (section IIII All to the realm of phase transition dynamics 
can be more expedient. Besides, nucleation does not commonly occur homogeneously in a matrix by means of solely 
thermal and concentration fiuctuations. It is induced most often by heterogeneities and defects in solids, e.g., grain and 
interphase boundaries, dislocations, stacking faults, free surfaces and vacancies or their clusters. These micro-domains 
facilitate or enhance the nucleation rate for formation of a new phase during quenching [23 |. 

IV. CONCLUSION 

The equations for nucleation, growth and coarsening of second phase precipitates are extended to account for 
the non-isothermal situations. The results of our simulation of the heat treatment of a Zircaloy-2 specimen clearly 
illustrates the infiuence of cooling rate during quenching on the properties of SPPs. A lower cooling rate results in 
larger precipitates with a smaller number density. This calculation is in qualitative agreement with the observations 
made on a similar kind of heat treatment of this material [^Jl, which also indicate that SPP size and density are 
quenching rate dependent and they impact macroscopic properties of the alloy. 
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FIG. 4: Evolution of Zr(Fe,Cr)2 SPPs in a Zircaloy-2 cylinder during the heat treatment. Diagrams, left to right, top to bottom, 
show calculated (i) mean precipitate radius, (ii) precipitate number density, (iii) nucleation rate, (iv) matrix supersaturation, 
(v) matrix solute concentration (Fe+Cr), and (vi) temperature in the first ten seconds. 
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APPENDIX A: NUCLEATION RATE CLOSE TO THE CRITICAL DROPLET RADIUS 



When the droplet (precipitate) radius R is close to its critical value Rc, the free energy, F(R) in Eq. ^ can be 
approximated by, F{R) = - 0.5(d - l)Sd<7R'^-^{SRf , with SR = R~ Rc'^. The integrand in Eq. Q is sharply 
peaked at Rc, thus the steepest descent method is used to evaluate the integral, resulting: 

J = {2n)-^/^Bc[(3\F"iRc)\y^\s{Rc), (Al) 

Substituting ns{Rc) — n/ exp [— /3i^(i?c)] and simplifying 

J = {2n)-^/^nceRcni exp {-[3Fc). (A2) 

Here, Vic = Bc(3\F"{Rc)\ is the nucleation frequency. Be = B{Rc) the critical kinetic coefficient for nucleation and 
e = {(3\F"{Rc)\Rl)~^''^ with F"{Rc) = [d'^F/dR^) ^^j^^. Utilizing Eq. CJ, we write, = SdBcl3(j{d - l)Ri-^ and 
— SdP<j{d — l)Rf~^ . Placing these last two relations into Eq. ljA2|l . we obtain Eq. © of the main text. 



APPENDIX B: STATISTICAL DESCRIPTION OF NUCLEATION RATE 



The starting point for the statistical mechanics approach to nucleation theory is the Ginsburg-Landau free energy 
functional for a scalar field variable ilj{f): 

(Bl) 

where ^ — ^{r) is the order parameter, e.g., for a binary alloy it is the local concentration, K \s a phenomenological 
constant designating the amplitude of the spatial gradient of -0 and U is the potential energy, which for a "t/;*" 
field theory is U{ili) = —^-0^ + ■^V'^ ~ ^"0- Here the coefficients tq and mq are, in general, functions of pressure 
and temperature and h is an acting external field, which for the case of binary alloys is equivalent to the chemical 
potential. The potential U has two minima, '4'±{h). For /i > 0, 0+ expresses the stable phase and ^- the metastable 
phase, whereas for h < the roles of 0+ and ip- are reversed (Fig. 1^. 



F(V.) 



dr 



(V0)2 + 
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To proceed, consider the grand-canonical partition function for the system, defined as: 

Z{h) = / S^Pe->^^'^^\ (B2) 



The extrema of the integrand in Z are the solutions of the equilibrium equation, i.e., 

if VV + roi} - UQip^ + h = 0. (B3) 

One solution of interest, for small h, is the saddle point solution, ip. It describes a radially symmetric particle of the 
nucleating phase embedded in a spatially uniform metastable background. Its profile is approximated by 



1 . , , s 1 



r-R 



-i^++^-) + -ii^+-ij-)tB.nh — =- , (B4) 



V2^ 



where ipif = 0) = ip- and ipif co) = ip+ are the stable and metastable values of the order parameter, corresponding 
to the two minima of U{ip), describing a stable particle in a metastable surrounding, r is the radial position, R the 
particle radius and ^ — {K/2roy^^ is the interface thickness called the correlation length. In a ^p^ field theory the 
correlation length is related to the capillary length £, cf. section III AL according to ^ = 6^ Q. 

The free energy F{^p) may be expanded around the saddle point 'ip. Using Dirac's concise notation, we write 

F{^/j) = F+^(u\M\u)+0{u^), (B5) 

where F — F{ip), u — ip — ip, and A4 ~ — isTV^ — C/"('0). We note that since "0 is a saddle point there is one negative 
eigenvalue which reflects the instability of the critical particle. Moreover, there are d zero eigenvalues associated with 
the translational modes of the particle in d-dimension. More expHcitly, as Langer [23| showed, the eigenvalues in the 
second order in h can be expressed by 



Ao/ 



{l-l){l + d-l) 
Rl 



2f2 



1 + ^(^)J' ^ = 0,1,2,... (B6) 



where Rc — {d — l)V2Kro / (Shy/uo) is the critical radius of the nucleus. Here, Aqo is the only negative eigenvalue, 
whereas Aqi is the c?-fold degenerate zero eigenvalue. The energy of the lowest mode ^ = is negative, implying that 
in the context of a steepest descent calculation the critical particle generates the imaginary part of the free energy in 
the metastable state. 

As it has been pointed out in in the case of interest {h 0, Rc oo), the eigenfunctions of M. form a 
band of soft modes which describe the surface excitations of the spherical critical particle. In fact, the eigenfunction 
corresponding to Z = 1 is the Goldstone mode associating with the spontaneous breaking of translational invariance 
by the center of the droplet. Apart from the unstable and translational modes the remaining eigenvalues of A4, A„; 
are positive describing stable distortion of the critical spherical particle [27| . 

To calculate the nucleation rate, we first express the partition function IIB2|l in the form of the Hnked-cluster 
expansion, viz.. 



Z = Zoexp(|i), (B7) 



where Zq represents the contribution from the metastable minimum and Zi that from the path ip. The contribution 
from the metastable minimum Zq = exp[— /3i^('0)], with ^/'(r) — tJj^ + ri[r) is written in the form 

Z„ = cxp[-/3F(v;+)] n ij^) "\ (B8) 

where A^*^^ are the eigenvalues of the equation: M.ii"f]f^ — ^f^vf^'' with Mq = — ifV^ — ro/2 + 1.5moV'+. Here 
■0+ = 2(ro/3uo)^^^ COS0 where 36 = arccos[1.5/i(3Mo/ri])^/^] is the the stable solution of Eq. ljB3ll for V^V = and 
h > 0. The saddle point contribution to the partition function can be written as 

' / ^ s(2/+l)/2 

Z,=exp[-m^)]JVn(^) , (B9) 

n.l 



11 



where the prime denotes that the zero modes are removed and J is the Jacobian that arises in integration and V 
denotes the contribution of zero modes and is proportional to the system volume V. The free energy of the system is 
T = limv/_oo (i^'^V^^ InZ. Writing Eq. ||B2|| in the form InZ = In^o + Zi/Zq and using Eqs. l|EHll and the 
droplet free energy is expressed by 



1^ 

T = exp[-/3F(V;)] J^Vr^/^ exp [ ^ In 



(BIO) 



where vi is the degeneracy of the c?-dimensional spherical harmonics, e.g., for d — 2, vi = 2 while for d = 3, i^i = 21 + 1. 
The upper limit of the sum L is given by L ~ R^ro, for which the approximation IIBfill becomes invalid. Giinther, 

Nicole and Wallace have calculated the factors in Eq. ljBlfl|l . Placing p = r^^^ / (|ft,|yuo), these factors are 



J4-d)/2 
'O 



Uo 



-p 



J ^ p 

and the spectrum of excitations are calculated to be 



_ (d-l)d/2 (i-d)d/i -d/2 



exp 



ip-'^exp(fcop^-i) d^i 
exp (^kop^^ d = 3 



ip 



-2/3 



(Bll) 
(B12) 

(B13) 



where c and are dimensionless constants and i = \/— 1. The nucleation rate a la Langer Q is J = (3k^{!F)/tt with 
K being a kinetic factor; Eqs. ijBllll . llB12p and ljB13|l can be assembled according to Eq. ljB10|l to yield the GNW 
relation [l3| in the regime of interest {h 0, Rc ^ oo), viz.. 



lim J = 



Ah 



At 



\ uo ) P 



exp 



^(■l-d)/2 



1/2 s 3/2 



3/2 
'^0 \~) 



p'^/^ exp 



1/2 



(cp'^-i + . . . ) d^3 
d = 3 



(B14) 



Now by transcribing the V'*-theory parameters to the pertinent variables in alloys, we identify cx^J^'^ p (xo/x), 

where ad = Tq* ''^^^/uo is an auxiliary variable and the variables x and xq are defined in section III Al Setting 

— 1/2 

£. = (with K — 2), we convert Eq. ljB14ll to 



d/(d 



d(d-3)/2 



7/3 



exp 



d^3 



(B15) 



From here on, to obtain Eq. ifTTijl is fairly straightforward, as it has been explicitly shown by Langer and Schwartz 



APPENDIX C: NUCLEATION DELAY TIME 



In order to calculate the time-lag for nucleation in classical theory, one needs to solve the full Fokker-Planck equation 
(PJ. Here, however, we can make a rough estimate of this quantity for the case of t >> B~^, meaning that; we are 
interested in time intervals much larger than the time of relaxation ((x B^^). We assume that n — n{R,t) obeys: 

(CD 

subject to the boundary conditions n{0,t) = m, n{s,t) = and the initial condition n(i?,0) = for s > Rc. The 
solution of this boundary value problem is expressed as 

n(R,t)^ni(l ) ^) — sin( )e ^ — . (C2) 

^ S ' TT ^-^ m \ s J 

m—l 
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The nucleation rate is related to the gradient of n{R,t), viz., 

BcUi 



i?4(0,g), (C3) 

_R=s 



where 'd4{u,q) = 1 + 2 X]m=i(— cos (2mu) is an elliptic theta function [2a|, here with q = exp (— i?c7r^t/s^) 
Employing the Courant-Hilbert identity (2fl| . 



Y: exp(-WC)(-l)" = ^ Y: exp[- "^"^~//'^ ], (C4) 



m— — OO m— — oo 



Eq. IIC3|l can be expressed in the form 

J(.,i) = ^(^)V2^,(0,e-^), (C5) 

■y/TTS t 

where td = /Be is the delay time. For short times, t << td, power series expansion of the theta function gives 

J{s,t) ^ ^^(y) ^ [e-- +0{e~ — )]. (C6) 

Equation lj("5p offers a zeroth approximation to the time-dependence of the nucleation rate, however, it illustrates the 
essential physics of the phenomenon. Detailed calculations of the time-dependent nucleation rate and the delay time, 
within the classical nucleation theory, can be found in Refs. j3nLl3ll|. 
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